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Abstract 

The Toda criterion of the Gaussian curvature is applied to calculate ana- 
lytically the transition energy from regular to chaotic motion in a schematic 
model describing the interaction between collective dipole and quadrupole 
modes in atomic nuclei. 
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In recent years great interest has been shown in nonlinear effects in 
nuclear 1 ) and particle physics 2 ). These phenomena may be considered the 
precursors of order-chaos transition. In these systems this transition has 
usually been studied numerically with Lyapunov exponents and Poincare 
sections 3 -*. Less attention has been paid to analytical criteria. 

In this paper we study the hamiltonian: 

H = l -{B iq \ + C iq l) + l -{B 2 q 2 2 + C 2 q 2 2 ) - k iq \q 2 , (1) 

introduced by Bohr and Mottelson 4 ) in order to describe the coupling between 
the giant-dipole resonance (GDR) and the quadrupole excitation, later used 
by the authors of reference 5 ). The vibration in the qi direction is associated 
with the GDR, while the q 2 motion plays the role of quadrupole deformation 
of the nuclear surface. Here the parameters of inertia Bi and of stiffness 
Ci (i = 1,2) are real constants. The perturbing potential k\q\q 2 induces 
a splitting of the dipole frequency caused by the static quadrupole defor- 
mation. Bohr and Mottelson 4 ) give the estimation k\ = (5/2)y / 3/(27r)Ci 
for the coupling parameter, in order to reproduce the experimental data of 
fotoabsorption in deformed nuclei 4 ). 

The aim of this work is to apply the Toda criterion of the Gaussian 
Curvature 6 ) to the Hamiltonian (1) in order to calculate the energy of the 
onset of chaotic motion, henceforth E c . 

The Toda criterion is based on a local estimation of the rate of separation 
of neighboring trajectories in the classical phase space of the model 6 )' 7 ). To 
obtain the time evolution of a dynamical system with the Hamiltonian 

2 



where pk = Bkqk, k = 1,2, the following equations have to be solved: 

rfq OH dp OH 



(3) 



dt <9p ' dt <9q 
where q = (51,52) an d p = (pi,P2)- The linearized equation of motion for 



the deviations are 



T(q(t)) 



where M^" 1 = 5ijB~ l , and 

where q(t) is the solution of (3). The stability of the dynamical system is 
then determined by the eigenvalues of the 4x4 stability matrix 

M- 1 \ 

. (6) 

-S(t) ) 

If at least one of the eigenvalues Aj of the stability matrix T is real, then the 
separation of the trajectories grows exponentially and the motion is unstable. 
Imaginary eigenvalues correspond to stable motion. 

To diagonalize the matrix T, we must first solve the equation of mo- 
tion (3). The problem can be significantly simplified by assuming that the 
time dependence can be eliminated by replacement of the time-dependent 
point q[t) of configuration space by a time-independent coordinate q, i.e. 
T(q(t)) = T(q). The eigenvalues then are 



A = ±[-b ± Vb 2 - 4c]^ , (7) 
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where 



1 2 ^ ~dqj ' 



d 2 V 8 2 V u » 



dql 
d 2 V 



(8) 

dql dql y dqiq 2 ' 1 ' ^ 
Now, if b > then with c > the eigenvalues are purely imaginary and the 
motion is stable, while with c < the pair of eigenvalues becomes real and 
this leads to exponential separation of neighboring trajectories, i.e. chaotic 
motion. The parameter c has the same sign as the Gaussian curvature Kq 
of the potential-energy surface: 

d 2 V d 2 V _ / d 2 V \2 
K G {qi, 92) = r; ; .oay.. , ,o2 V .„2 ■ (10) 



[i + (^) 2 + (Sf) 2 f 

For our nuclear model the potential energy is given by 



V(qi, q2) = 7>( Ciq2 + ° 2q ^ ~ kiq2q2 



and the T matrix reads 



(11) 



r(gi,92) 



V 

























Ci + 2fcig 2 










2fci5i 


-c 2 









/ 



(12) 



The potential has one minimum for q 1 



Q2 



with V = 0, and two 



saddle points for q 1 = ±y^^, g 2 = These are the points for which 

the curvature criterion is exact (they are the fixed points of the Hamilton 
equations). The origin is a stable elliptic point and the two saddle points are 
unstable hyperbolic points. 



The dissociation energy of this model is Ed = -^r- and it occurs at the 
saddle points of the potential energy surface. Classically, a particle located 
within the potential well with energy < E < E D is always bound. 

The Gaussian curvature vanishes at the points that satisfy the equation 

C 2 (C 1 -2k 1 q 2 )-4klql = 0, (13) 

where 

9? = ||(Ci - 2*ift) , (14) 

with q 2 < At low energies, the motion near the minimum of the potential 
energy, where the curvature is positive, is periodic or quasi-periodic and 
is separated from the region of instability by a line of zero curvature; if 
the energy increases, the system will be, for some initial conditions, in a 
region of negative curvature where the motion is chaotic. In accordance 
with this scenario, the energy of order— ^chaos transition E c is equal to the 
minimum value of the line of zero Gaussian curvature K G of the potential- 
energy surface of the system 6 )' 7 ). The energy of the zero-curvature line is 
determined by the expression 

V(K G = 0, g2 ) = C 2q i-^ q2 + ^. (15) 

It is easy to show that the minimal energy on the zero-curvature line is given 

by 

£c = V(^ = 0, 92 = ^) = ||. (16) 

This is the energy of the onset of chaos of the model and we see that E c = 
E D /2. 
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With the estimation k\ = (5/2)y2/(27r)Ci we have 

By using the charged liquid drop approximation 4 ) for the stiffness coefficient 
C 2 , we obtain for E c the following relation 

^ ,b s .2 3 Z 2 e 2 ,71 , . 

^ = <7- 4i -lte^4W»150' < 18 » 

where b s ~ 17 MeV, r ~ 1.2 fm, A and Z are the atomic and proton 
numbers respectively (see Figure 1). With the liquid-drop approximation 4 ), 
the Coulomb repulsion counteracts the effect of the surface tension and leads 
to negative values of C 2 for Z 2 /A > (10/3)(6 s r /e 2 ) ~ 49. 

We can use the experimental quadrupole frequency hu 2 and B(E 2 ) values 
to determine 4 ^ the restoring force parameter C 2 . In this case 

E c = ^hcu 2 (^-ZeR 2 ) 2 B(E2; ^ 2)" 1 . (19) 

The experimental values of C 2 exhibit variations greater by a factor 10 above 
and below the liquid-drop estimation. This effect can be understood in terms 
of the approach to instability as particles are added to the closed shells and 
so E c is in the range ~ 0.5-5 MeV depending on the nucleus. 

The curvature criterion is able to characterize the local behaviour of the 
system (for example the local instability) and may give only a signature of 
the global properties (e.g. the global instability). As is well known, a very 
useful tool for the study of global properties is provided by the Poincare 
sections 9 ^. With this aim the classical trajectories have been calculated by 
the fourth order Runge-Kutta method. In order to avoid numerical errors 
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connected to the use of finite temporal intervals, a first-order interpolation 
has been used 10 ) . 

The Hamilton equations of the systems are: 

qi = Bipi , q 2 = B 2 p 2 , 

Pi = -Ci5i + 2hqiq2 , P2 = -2C 2 g 2 + hql • (20) 

Figure 2 shows the Poincare sections for different values of the energy. Be- 
low the critical energy E c the system is regular but above E c , the Poincare 
sections clearly show an order-chaos transition as the energy increases. 

In conclusion, the Poincare sections show that in our schematic model, 
which describes the interaction between collective dipole and quadrupole 
modes in atomic nuclei, the Toda criterion of the Gaussian curvature gives 
a quite good estimate of the energy E c of the onset of chaos. This critical 
energy E c depends on the stiffness coefficient of the quadrupole vibration, 
which can be estimated analytically within the liquid-drop approximation 
or can be experimentally determined by using the quadrupole frequency hu 2 
and the B(E 2 ) transition probabilities. 

LS is greatly indebted to Prof. M. Robnik and Dr. A. Rapisarda for 
many enlightening discussions. 
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Figure Captions 



Figure 1: The critical energy E c versus the atomic number A, with Z = A/2. 
Figure 2: The Poincare sections for different values of the energy; from 
the top-left, clockwise: E = 0.02, E = 0.03, E = 0.035 and E = 0.04. 
We chose B 1 — B 2 — 1, C\ — 100 and C 2 = 1. The critical energy is 
E c = ^C 2 ~ 0.021. 
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